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Abstract. The recent observations of type la supernova strongly support that the 
universe is accelerating now and decelerated in the recent past. This may be the 
evidence of the breakdown of the standard Friedmann equation. We consider a general 
modified Friedmann equation. Three different models are analyzed in detail. The 
current supernova data and the Wilkinson microwave anisotropy probe data are used 
to constrain these models. A detailed analysis of the transition from the deceleration 
phase to the acceleration phase is also performed. 
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1. Introduction 

The recent observations of Type la supernova indicate that the expansion of the universe 
is speeding up [T]. Observational resuhs also provide the evidence of a decelerated 
universe in the recent past On the other hand, the cosmic background microwave 
(CMB) observations indicate that the universe is spatially flat as predicted by the 
inflationary models A dark energy component with negative pressure behavior 
which dominates the universe, was proposed to explain a flat and accelerating universe. 
One simple candidate of dark energy is the cosmo logical constant. However, there are 
some problems with the cosmological constant although cold dark matter cosmological 
constant models are consistent with the current observations. Why is the cosmological 
constant so small and not zero? Why does the cosmological constant become significant 
now? The quintessence models avoid some of the problems There are also other 
models, like tachyon filed as dark energy [5 . But the property of dark energy is still 
mysterious. One logical possibility is that the standard Friedmann equation may need to 
be modified. In this scenario, the universe is dominated by ordinary pressureless matter, 
but the law of gravity and the standard Friedmann equation are modified. The idea of 
modifying the law of gravity is not new. The modified Newtonian Dynamics (MOND) 
was first used to explain the rotation curve in place of dark matter [7]. In MOND, 
the Newtonian gravity M/r^ is replaced with M/r^ + y/M/r. Since M/r'^ gives the 
standard Freidmann equation ~ p, M/r'^+y/M/r may provide a modified Friedmann 
equation ~ p^/^ In p + p^^^ j2j . Recall that the brane cosmology gives a non-standard 
Friedmann equation ~ P+P^ [HHHOni- Along this line of reasoning, Freese and Lewis 
recently proposed the Cardassian expansion in which the universe is dominated by the 
ordinary matter and the Friedmann equation becomes ^ p + p^ jTTj. The Cardassian 
model was later generalized to a more general form ~ g{p) [12]. In addition, several 
authors modified the Friedmann equation as if ^ + ~ p motivated by theories with 
extra dimensions ^3 In this paper, we first consider a model which is equivalent 
to the generalized Chaplygin gas model ^Hj in terms of dynamical evolution. Then 
we consider the generalized Cardassian model. At last we consider a model proposed 
by Dvali, Gabadadze and Porrati (DGP) ^Hl- We use the ten new supernovae at 
z = 0.36 — 0.86 jini and the the first year Wilkinson Microwave Anisotropy Probe 
(WMAP) temperature (TT) and temperature polarization cross correlation (TE) data 
[T7] to constrain these models. We also investigate the transition from the decelerated 
phase to the accelerated phase. 

For a spatially fiat, isotropic and homogeneous universe with both an ordinary 
pressureless dust matter and a minimally coupled scalar field Q sources, the Friedmann 
equations are 




(1) 



a 



(Pm + PQ + Spq), 



(2) 



a 
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Pq + 3H{pq+pq)=0, (3) 

where dot means derivative with respect to time, pm = Pmol'^o/'^)^ is the matter energy 
density, a subscript means the value of the variable at present time, pq = Q'^ /2 + V{Q), 
Pq = — V{Q) and V{Q) is the potential of the quintessence field. The modified 

Friedmann equations (MFE) for a spatially flat universe are 

H' = Hlg{xl (4) 



a 



3HqX , . f p + p^ 



H',9{x)---^g'{x)^^j, (5) 

p + 3H{p + p) = 0, (6) 

where x = SvrGp/Sifg = xo{l + zY during the matter dominated epoch, 1 + z = a^/a is 
the redshift parameter, g{x) = x + ■ ■ ■ is a general function of x and g'{x) = dg{x)/dx. 
From phenomenological point of view, four dimensional gravity is modified so that we 
get a general function g{x). Chung and Freese argued that a general g{x) is possible 
if our universe as a three brane is embedded in five dimensional spacetime For 
example, g{x) ~ x + in Brane cosmology, g{x) ~ x + x" in Cardassian model. On 
the other hand, we can think the additional terms g{x) — x as dark energy component. 
For the Cardassian model, the additional term x" can be mapped to a dark energy 
component with constant equation of state parameter ujq = n — 1. In general, we get 
the following relationship between the dark energy equation of state parameter and g{x) 
xg'{x) - g{x) 
g[x) - X 

The equivalent dark energy potential can be found from the following equations 

= p^[g\x) - II (8) 

y{Q) = - 0.5x - 0.5x(7'(x)]. (9) 

For instance, if g{x) ~ x + x", we find that V{Q) ~ [sinh(AQ + Note 
that the universe did not start to accelerate when the other terms in g[x) started to 
dominate. The linear density perturbation of this model is given by 

5 + 2H5 = A7rGp6[g{x) + 3x/(x)]. (10) 

For the matter dominated fiat universe, p = Pm and p = Pm = 0. Let i^mo = 87iGpo/3Hq, 
then xo = ^mo, g{xo) = 1. In general, x = ^^mo(l + + firo(l + z)^, where 
Qj-o = 8.35 X 10~^ is the current radiation component |18j . 

2. Analytical Method 

The location of the m-th peak of the CMB power spectrum is parameterized as |I9| 

lm = {m- (pm^A, (11) 

where the acoustic scale Ia is 

, vr To - _ TT /cf-((7Ko(l + + ao(l + z)'])-'/'dz 

^ c- ris Csmg[n^M + zr + n,o{l + zr])-y^dz' ^ ^ 
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Cs = 0.52, the conformal time at last scattering and today tq are 

ri^ , f°° dz , , 

ris= dT= , (13) 



ro - r dr ^ r (14) 

^0 ^0 aoHo^g[n^o{l + zf + n,o{l + zy] 

The recent WMAP results give the positions of the first two acoustic peaks as Ipi = 
220.1 ± 0.8 and lp2 = 546 ± 10, respectively [T7j. The third peak is given by the 
BOOMERanG measurements as Ips = SAbtH [SO!- We first use equations (HU-jlll) with 
01 = 0.3074, 02 = 0.2819 and 03 = 0.341 to constrain the models considered below, 
then we use the full 1350 WMAP TT and TE data [HI by a modified CMBFAST code 
version 4.5.1 j2Tj to constrain the parameters. For the fit to full WMAP data, a scalar 
power spectrum with normalization 0.833, spectral index 0.93 and running index slope 
-0.031 are assumed. Other cosmological parameters are chosen as follows: h = 0.71, 
= 0.044, Tcmb = 2.725, Helium abundance Ine = 0.24, Number of massless neutrinos 
is 3.04 and g* = 10.75. 

The luminosity distance di^ is defined as 

, . . . , r^o dt' 1 + z du , ^ 

d^{z) = ao(l + z) -— = —— / , (15) 



The apparent magnitude redshift relation becomes 

m{z) = M + 5 logio M^) +'25 = M + 5 logio '^l{z 

du 



M + 5\og 



10 



(1 + ^) r 

JO 



(16) 



'g[^lmoil + uY 

where T>i^{z) = Hodi^{z), M is the absolute magnitude and A4 = M — 51ogj^Qifo + 25. 
The nuisance parameter A4 can be determined from the low redshift limit at where 
T>i^{z) = z. The parameters in the models are determined by minimizing 

i 

where cxj is the total uncertainty in the observation. We use equations flT^ and (fT^ to 
fit the ten new supernova data with host galaxy correction ^H] • The nuisance parameter 
M. is marginalized when we fit the supernova data. 

The transition from deceleration to acceleration happens when the deceleration 
parameter q = —d/aH"^ = 0. From equations and we have 

3 

^[^1^0(1 + Zg=o)^] = -l^mo(l + 2;g=o)V[^mo(l + ^g=o)^], (18) 

3 

go = 2^m05''(^mo) - 1- (19) 



Constraints on alternative models to dark energy 



5 



3. Chaplygin Gas Model 

In the framework of MFE, the generahzed Chaplygin gas model pc = —A/p'^ becomes 

g{x) = x + nQo[As + (1 - As){x/n^of]'^^, 

where figo = 1 - ^mo - ^ro, f3 = I + a and As = {^-kG /^HIVLqqY A. To recover the 
standard Friedmann equation at early times, we need Ag ~ 1. When A^ = 1, the model 
becomes a standard A-model. Now we have 
du 



du , . 

^f{u) + ^^Qo{^ + (1 - ^)(i + uf + ao(i + uY/n^oYY/p 

where /(m) = VI^qY + + ^ro(l + uY- Since g\x) = 1 + VIqqY - As)[As + (1 - 
As){x/VL^QYY''^''^{x/VLri,QY ^ together with equations (HHI) and (HH)), we have 

^Y + ^,=oY[As + Y- AYY + z,=oY'Y-"' = ^ - ^(i - + ^,=o)'^ (21) 

go = ^-^A,(l-r]^o), (22) 

The first three peaks in CMB power spectrum favor a cosmological constant model 
with As = I 01 (3 = I. The best fit to the WMAP TT and TE data is Vl^q = 0.26, 
As = 0.999 and (3 = 1.43 with = 1448.3. The ten supernova data also favor a 
cosmological constant model with As = 1 or (3 = 1. By using the best fit parameters 
to WMAP data, we get Zg=o = 0.78 and go = —0.61. Because As = 0.999, so g{x) ^ x 
at early times. The best fit result to the supernova data is shown in figure ^ and the 
WMAP TT power spectrum with the best fit parameters is plotted in figure El 



4. Cardassian Model 

We take the generalized Cardassian model 

g{x) = x[l + 

where B = (fi^g ~l)/^mo"~^\ a > and n < 1 — 1/3(1 — J^J^q). At early times, g{x) ~ x, 
so the standard cosmology is recovered. When n = 0, gY) = B^^^'Y + x°' / BY^°' ■ For 
the special case a = 1 and n = g{x) = x + B which is the standard cosmology with a 
cosmological constant. If we take a = 1 and n = 1/2, then we have gY) = x + B^x. If 
we think the generalized Cardassian model as ordinary Freidmann universe composed 
of matter and dark energy, we can identify the following relationship for the parameters 
in the Cardassian and quintessence models 

(n-YY-^'^o) 
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Figure 1. The Best fit curve to the supernova data. GeneraUzed Chaplygin gas: 
^^mo = 0.27, As — 0.94 and (3 = 1.0005; Generalized Cardassian model: f7niO ~ 0, 
a = 0.23 and n = 0.49; DGP model: n^o = 0.21 



The generalized Cardassian model gives 



Jo 



du 



g[n^Q{l + uf + + u) 



du 



(23) 

(24) 

(25) 
(26) 

The best fit parameters to the first three peaks in CMB power spectrum are: fimo = 0.28, 
n = —0.032 and a = 0.5 for the generalized Cardassian model and f2mo = 0.26 and 
n = —0.12 for the Cardassian model. For the Cardassian model a = 1, the best 
fit prameters to the full WMAP TT and TE data are Q^o = 0.25 and n = 0.02 with 
= 1461.0. The best fit parameters to the ten supernovae data are: flrno ~ 0, n = 0.51 
and = 6.66. For the generalized Cardassian model, the best fit parameters to the 
ten supernovae data are: flmo ~ 0, n = 0.51, a = 0.54 and = 6.66. In fact, the data 



/o {f{u)[i + {Q-^ - i)[(i + uf + ao(i + w)v^]^o]"("-^)]^/"}'/' ■ 

Combining equation (j23|) with equations (fT8|l and (fT9|) . we get 

l + z,=o = [(l^J-l)(2-3n)]V3"(i-), 



1 3, 
- + -(n 

2 2^ 



Constraints on alternative models to dark energy 



7 



is not sensitive to n and a. We get almost the same value of with different a and 
n. So the supernova sample is too sparse to determine the parameters in generalized 
Cardassian model. By using fl^o = 0.25 and n = 0.02, we get go = —0.6, Zq=o = 0.82 
and uqq = —0.98. The best fit curve to the supernova data is shown in figure [T] The 
CMB TT power spectrum with the best fit Cardassian parameters is drawn in figure El 



5. DGP Model 



The model is g{x) = [a + \/a? + xY |T3^, where a = (1 — Q^o)/2. At high redshift, 
g{x) ~ X. So the standard cosmology is also recovered in the early times. The conformal 
time is 

du r du ^27^ 







For this model, we find that go and the transition redshift Zq=Q from decelerated 
expansion to accelerated expansion are 



1 + Zq=Q 



2(1 -VL 



mOj 



\2 



mO 



1/3 



(2^ 



= TTT^- ^^^^ 

In this model, we have only one free parameter f2mo- The locations of the first three 
peaks in CMB power spectrum give that fimo = 0.179. The best fit parameter to the 
WMAP TT and TE data is fi^o = 0.174 with = 1485.6. The best fit parameter 
to the supernova data is Q^no = [0, 1.0] centered at 0.21 with = 6.67. If we map 
the additional term in the right hand side of the Friedmann equation to dark energy 
component, we have a;Qo = —1/(1 + ^mo) = —0.83, go = —0.48 and Zg=o = 0.81 if we 
take Qrao = 0.21. The best fit curve to the supernova data is shown in figure C] and the 
best fit curve to the WMAP data is shown in figure |21 



6. Discussions 



One attractive feature of the generalized Chalygin gas model is that it can be considered 
as a unified dark matter and dark energy model. In the unified scenario, both the dark 
energy and dark matter components are modelled by the generalized Chaplygin gas. In 
this paper, we are interested in modelling the generalized Chaplygin gas as dark energy 
only. We found the best fit parameters in the generalized Chaplygin gas model are 
^mo = 0.26, As = 0.999 and f3 = 1.43. Cunha, Alcaniz and Lima found that > 0.73 
by using the combined x ray data of galaxy clusters and supernova data and taking a 
prior fiino = 0.3 ^Hl- Amendola et al found that 1 < (3 < 1.2 and 0.8 < < 1 by using 
WMAP TT data In the context of MFE, the generalized Chaplygin gas model 

tends to be the A-model in order to recover the standard cosmology at early times. 
Our results show that the generalized Chaplygin gas model is almost the same as the 
A-model. By using the best fit parameter, we get Zq=o = 0.78. 
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Figure 2. The Best fit curve to wmap data. The dashed Hne is for generalized 
Chaplygin gas model with fimo = 0.26, Ag = 0.999 and (3 = 1.43. The Solid line is for 
Cardassian model with fimo — 0.25 and n — 0.02; The dash dotted line is for the DGP 
model with fimo = 0.174. 

For the Cardassian model, the WMAP TT and TE data give that Q^o = 0.25 
and n = 0.02. The 10 supernova data gives that fl^no ~ and n = 0.51. The error 
from supernova data fit is also large. Zhu and Fujimoto also found low Q^^q from 
supernova data 22 . Sen and Sen found that 0.31 < n < 0.44 and 0.13 < Q^no < 0.23 
from supernova data and peaks in CMB power spectrum, n < 0.3 from WMAP data 
[23] • Godlowski and Szydlowski found that n^o = 0.48lo:?3 and n = -OAtiZ, and 
^mo = 0.5llo;og and n = —1.2lii'])l from different sets of supernova data jSlj. Frith 
found that 0.19 < Qrao < 0.26 and 0.01 < n < 0.24 from supernova data [23]. In 
general, different analysis gives different results, but they are still consistent with each 
other at 99% confidence level. Dev, Alcaniz and Jain found that the best fit to the 
gravitational lensing effect is n = 0.76 and a = 2.4 by assuming Q^q = 0.3 for the 
generalized Cardassian model j^. Taking the best fit result to WMAP data, we get 
Zg=o = 0.82. 

For the DGP model, the best fit parameter to the WMAP TT and TE data is 
^mo = 0.174 and the best fit parameter to the ten supernova data is fimo = 0.21. 
Deffayet et al found that f2mo = 0.18to;og from supernova observations and f2mo = 0.3 
from CMB constraints fI7\. However, the constraint from supernova data and CMB data 
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are still consistent with each other at la level. If we take the bigger value flj^o = 0.21, 
then we get Zq=o = 0.81. Multamaki, Gaztanaga and Manera considered the effects of 
Cardassian model and DGP model on large scale structure and found that these effects 
are different from that of A-model |2H]- 

Due to the uncertainties in the observational data, it is still difficult to discriminate 
different dark energy models and alternative models. Wang et al showed that future 
supernova data such as SNAP | could differentiate various models if Q^no is known with 
10% accuracy [SH]- On the other hand, if we can determine the dynamical evolution of 
equation of state parameter of dark energy with high precision using model independent 
method, then we can discriminate different models. Alam et al showed that future SNAP 
data will provide important insights on the nature of dark energy at high redshifts |30j . 
Turner and Riess showed that Zq=o ~ 0.5 Daly and Djorgovski found that Zq=Q > 0.3 
|31j . In conclusion, the above models are consistent with current observations. 

In terms of model building, we can assume a particular scale factor which manifests 
early deceleration and later acceleration, then find out the form of g{x). For example, 
we take a{t) = ao[sinh(t/to)]'^/a, then we find g{x) = Ax^^^^ + B. This model recovers 
the standard cosmology when (3 = 2/3. 
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